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Abstract 

The equilibration of two coupled reservoirs is studied using a Green func- 
tion approach which is suitable for future development with the closed time 
path method. The problem is solved in two parameterizations, in order to 
demonstrate the non-trivial issues of parameterization in both the interme- 
diate steps and the interpretation of physical quantities. We use a covariant 
approach to find self-consistent solutions for the statistical distributions as 
functions of time. We show that by formally introducing covariant connec- 
tions, one can rescale a slowly varying non-equilibrium theory so that it ap- 
pears to be an equilibrium one, for the purposes of calculation. We emphasize 
the importance of properly tracking variable redefinitions in order to correctly 
interpret physical quantities. 



1 Introduction 

Effective field theory is an increasingly important tool in the repertoire of theoretical 
physics. The problem of the equilibration of statistical systems has been studied 
in many contexts |2[ |[ ^, H, || , but has not been formulated in the language of 
effective quantum field theories, where re-parameterizations of basic quantities can 
obscure the physics of the process. 

Suppose we consider two systems with different mechanical and thermodynam- 
ical properties and couple these together. What is the outcome? Intuition tells us 
that such a system must achieve a state of thermodynamical equilibrium. However, 
it is not clear how this equilibrium is achieved using the language of effective field 
theory. In terms of particles, it is straightforward to picture the exchange of energy 
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between heat baths: there is a physical flow of point-Uke objects carrying momen- 
tum, which are exchanged until, eventually, there is an equal spectrum of energies 
on both sides of the interface. In a field theory however, one deals with quanta or 
modes of the field, so there is no direct notion of a flow of hot into cold or vice 
versa, only the occupation of states in a Fock space. 

The procedure of identifying a functional effective field theory involves many 
variable redefinitions and partial summations. Each re-parameterization distances 
the system from any intuition one might have about its behaviour. Our goal in this 
paper is to understand the process of equilibration in the language of field theory 
in the context of a specific model. We consider a pair of essentially free scalar 
fields with different masses and different initial temperatures which are coupled 
minimally at an initial time. For definiteness we assume that the systems are in 
direct contact, rather than interacting through a boundary layer, so that the two 
systems essentially occupy the same space. Although it is convenient to focus on a 
specific model, our aim is to extract general principles about non-equilibrium field 
theory. We employ the methods of general covariance during field redefinitions, in 
order to track the meaning of the quantities which arise in the formulation, and 
also to find the optimal parameterization for solving the equations. As is often the 
case with covariance, this leads to important physical insights. 

The notation and conventions are those of reference ^, with n spatial di- 
mensions and one time dimension. To simplify notation we shall use Schwinger's 
shorthand for the measures: 



/ 



(27r)"+i 



{dk) 
(dk) 



(27r)" 

dV^ = dtd'^x = (dx) (1) 

and so on. 

A non-equilibrium system remembers its history; it is not translationally invari- 
ant in time. It is convenient to parameterize Green functions in terms of even and 
odd variables over any interval from x to x': 

X — [x ~ x') 

X = + (2) 

The odd variables x characterize translational degrees of freedom, while the even 
variables x represent inhomogeneous impediments to translational freedom, or alter- 
natively, characterize local curvature in the dynamical and kinematical parameters. 

The formulation we adopt is based on [H, which is built around Schwinger's 
closed time path quantum action principle. We do not need to formulate the full 
closed time path problem here, but our approach is easily generalizable for future 
problems. We concentrate on understanding the basic unitary processes which act 
in the system, and employ the language of sources and Green functions to describe 
the development. 



2 The model 

Consider two systems which are brought into contact, so as to interact through a 
minimal coupling which is not manifestly Z2 {4> —(f) invariant. Once coupled, 
the systems influence one another. Since the systems are statistical mixtures of field 
modes, this influence can be viewed at two levels: as the microscopic interaction of 
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individual field solutions between the two systems, or as a change in the statistical 
averages of the two systems as a result of the microscopic interaction. It is at this 
latter level that one expects to achieve equilibrium. The model we choose begins 
as a free field theory, composed of a thermal mixture of quanta. This initial state 
is clearly idealized since one normally requires some form of interaction in order 
to achieve thermal equilibrium in the first place. In the interests of simplicity, we 
suppress all such non-linear interactions and assume that the initial state is two 
effectively free thermal baths. 

We formulate the physical system using the Schwinger action principle. The 
uncoupled systems have the action operators: 

51 = J{dx)[{d^^ind^^i) + ml} 

52 = j{dx){{d^^2nd^^2) + ml}. (3) 
Before coupling the equations of motion give us the eigenvalue equation: 

^^ = - f "°rt"' V tO=0- (4) 

V -D +mi J \ <i>2 J ^ ^ 

The Green function Q for the operator O satisfies the equation OG = I. We 
write the Green function as the sum of a particular integral and a complementary 
function 

g = G + (5) 

which satisfy 

dG = I; dT = Q. (6) 

In most pure state cases one sets T = Q and chooses a special solution for the 
Green function G which satisfies certain specific boundary conditions (retarded, 
advanced, etc.). The complementary function T becomes important in statistical 
many particle systems where mixtures of fields are significant. However, it is not 
possible to associate specific boundary conditions with the free particle solutions 
that are contained in the complementary piece (free particle solutions do not respect 
time-specific boundary conditions, since plane- wave solutions have no beginning or 
end) . The coefficients of plane wave solutions must therefore change in time in order 
to follow the development of the system, in a non-equilibrium many particle system. 
This time dependence is discussed in detail in the next sections. We obtain, 

G= \ -^r^ ? \ (7) 




and 



fi{k^)5{xi) \ 

V /2(fco)5(X2) )' ^ ^ 

Although nothing formally precludes off diagonal terms in where the upper right 
contribution is proportional to 6{xi) and the lower left contribution is proportional 
to 6{x2), the boundary conditions tell us that there is no interaction between the 
two systems, so off-diagonal terms would not make sense. 
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Once coupled, the fields are influenced by their mutual interaction. Since this 
interaction leads to modifications to the constraints (mass shell) of the systems, we 
denote the coupled variables without carets: 

51 = J{dx){{d^,<tiy{d''^i) + ml^i<^i + j*<^i + Ji^i} 

52 = J{dx){{d^^2nd''^2)+rhl^*2^2 + J2^2 + J2^2} (9) 

Both systems interact only through linear perturbations. In order to couple the 
systems together we choose: 

Ji = ml $2 

J2 = ml^i, (10) 

where nip is a squared polaron mass introduced for dimensional consistency so that 
each field is a source for the other. Starting from the total action, 

Stot = {Si + S2)\j^^^,^^j^^^,^^, (11) 

the operator equations of motion are, 

(-□ +mf)$i = -ml^2 

(-□+m^)$2 = -ml^i. (12) 

Note that the coupling destroys the diagonal structure of the equations of motion. 
We have: 

0$ = J$ (13) 

With O = O — J we can write, 

O$ = or f-°+5 ~1^^2)(t)=0. (14) 

In order to find the corresponding Green function we diagonalize the matrix O. 
The technique is standard. Wc solve the eigenvalue equation 

-° +2^' -"1^ 2 V ^0 = (-□ + ml) (^'). (15) 

and compute the invertible unitary matrix of eigenvectors U. The diagonalized 
matrix is obtained from the similarity transformation, 

O' = (C/1')-iC'[/-i = C/OCZ-i (16) 

The Green function matrix transforms according to the same rule as the operator 
itself: 

g ^g' = ugu-\ (17) 

which preserves the relationship 

Og = O'g' = I (18) 
The similarity transformation corresponds to a field transformation of the form, 

= UaB ifB, (19) 
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where A,B^1, 2. 

We proceed as follows. Solving the eigenvalue equation (|l5|) for the eigenvalues 
gives, 



= ± y rh^ + rrip 

— ™2 _i_ 2 



±m% . (20) 



with 



:= i(mi +TO2) ; ifi^ ^{mj - ml) (21) 
Using the identities, 

r 2 ~2\ ~ 2 I 2 

(m^ — rri]^) = —m ±171^ 

(m± - 7712) = TO^ ± toIj , (22) 
the invertible unitary matrix of eigenvectors is easily calculated: 

IT _ ^ ( "^P ("ifl - "1^) y ^ cos6l sin^ \ /23) 

/m* + (m2 _ ^2^)2 I -{ml-m') ml J' {-sine cosOj^ 

Using (|l|), dH) and (H) we obtain the diagonalized form of the operator O: 

0'=( -°+'^+ ,2 V (24) 

Note that when toi — 7112 (or m = 0) this result reduces to a familiar 7r/4 rotation 
of field variables. The matrix (23) becomes. 



which we use in (19) to obtain the simple rotational form of the transformation, 

$2 = ^(<^2-^l) (26) 

To rewrite the Green function in the eigenbasis of the coupled system we solve 
the equation O'Q' = /. As before, we separate the Green fimction into the sum of 
a particular solution and a complementary function: Q' = G' + 2-KiT' which satisfy 

O'G' = / ; O'T' = (27) 

The inverse of the diagonal operator O' is trivial to evaluate. We obtain 



1 



G' = ( ] = -°r^+ 1 I . (28) 







^21 ^22 / \ 



-□- 



In the coupled system, the diagonal form of the complementary function does not 
have to be preserved. However, we know that: 

T' = UTU-\ (29) 



5 



In addition, we know that T' must satisfy ( |27| ) which means that it must have the 
form, 

_( fiiiko)S{x+) fi2iko)S{x+) \ (^^. 

We assume that the original J- is diagonal just prior to the interaction begin turned 
on, so that if the interaction is turned on suddenly, remains diagonal: 

^ //r(Wxr) \ 

V f2{ko)S{x2) J 



where Xi and X2 are as yet undetermined mass shells. Using (|2^) and ( plD in (g9|) 
and comparing with (^0|) we obtain the relationship between the components of 
and JF': 

fll /l 2 /) , /2 • 2 /J 

= — cos -\ sm 

UJ^ UJi UJ2 

Ik = A cos^ + A sin^ 

fl2 / /2 /l \ n ■ n 

— ( jcost'smt^ 

UJ^ UJ2 ^1 

(A _ A)cos6'sin6'. (32) 

CiJ- U}2 OJl 



3 Covariant analysis at t 7^ 

Now we consider the behaviour of the system for i 7^ 0. We make use of the covariant 
method of [|| to solve for the Wightman or complementary Green functions for the 
effective theory. The off-diagonal terms are irrelevant to further analysis; they 
trivially satisfy the equations of motion by virtue of the fact that their particular 
solutions are identically zero and their complementary parts contain delta-function 
constraints. For notational simplicity we therefore refer to the diagonal components 
as G'li G;, G'22 G'_, f{, /; and /^^ := /-■ 

Our strategy is to find the dispersion relation for the positive frequency Wight- 
man function. This dispersion relation does not have a unique form: it depends on 
the ansatz used for the Green function itself. By using a conformally covariant form 
for the ansatz and introducing vector fields A^, we obtain a form for the dispersion 
relation which is independent of x. This approach is useful in two ways: it results 
in a neater form for the dispersion relation, and it allows the internal conformal 
symmetry to be seen directly and used to simplify calculations where the dispersion 
relation is not x independent. Note that we are not claiming that the theory is 
conformally invariant. We are only noting that an overall x scaling of the action 
can be used to express any statistical theory in terms of dimensionless variables, 
provided partial derivatives are replaced by covariant ones |8|. By formally intro- 
ducing covariant connections, we can exploit our ability to rescale a slowly varying 
non-equilibrium theory so that it appears to be an equilibrium one, for the purposes 
of calculation. 

The positive frequency Wightman functions satisfy the equations of motion: 

(-□ +ml)Gi+y{x,x') = 0, (33) 

where the index a takes values {+,—}. We attempt to find solutions to these 
equations using the ansatz, 

Gi+y{x,x') = 27rz /(dk)e*^-^+'r^-^^''iiA^%MlI». (34) 
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The formal dependence of the distribution functions and on the average coor- 
dinate X seems pecuhar in a free theory, where there is no apparent cause for a 
re-distribution of quanta amongst different modes in either time or space. Such a 
change would have to be represented by an additional source or sink (a potential), 
which we do not have. However, there is nothing in the field equations to exclude a 
combination of solutions which oscillates with a steady periodic signature. Such a 
steady-state oscillation could be fed in as a boundary condition. Another example 
would be the analysis of low-order terms in an interacting theory, where the effec- 
tive mass was changing with time m^it). Such a change makes no sense without a 
broader theory to explain it, but as an effective model it can still be analyzed in a 
consistent framework. 

The purpose of the auxiliary field A^^ is as follows. Since external boundary 
conditions and renormalizations might impose a time-dependence in the disper- 
sion relation, one needs a generic term to absorb the effect of the changing 
distribution. Its introduction allows us to formulate the theory in terms of a time- 
independent frequency, by moving the dependence of the average coordinate into 
the phase. 

Taking the derivative we obtain. 



d,Gi+yix,x') 



2ni I (dk) ik^, + -Aa,, + i I {df,Aa^)dX'' 



Gi+y{k,: 



(35) 



where, 



df,uja{x)/uja{x); F' 



1 + /^ 



(36) 



and 



X 

□ 



□ + d^d" + in 



(37) 



Notice that the second term of the first line above behaves as an auxiliary field. It 
is this term that will be absorbed by the auxiliary field. 

We proceed as follows. We make the choice Aafj = id^uja- Substituting into 
i Jq {di_iAav)dx"^ and integrating by parts we obtain, 



* / id^Aa.)dx'-' = i5"[uJa{x) - ^a(O)]. 



(38) 



When we substitute this result back into the dispersion relation we obtain. 



d-^Gi+y = 2m / {dk) 



^(ka^ + SJ'[u;a{x)-u:am) + liK, 



Gi+y (39) 



We re-express this result covariantly by defining the vector fca^(O) = i^a(0)(5^o+fc*^ip 
which allows us to write 



d^Gi+y = 27n J (dk) [^k,^{0) + ^{F^^^ - f^.^)]Gi+)' 



(40) 



Performing another derivative and substituting into (p3) we obtain. 



1 



ka.mm - ml + -[B^^iF!,^ - Qa,.) + [F', - n^f] + ik^iKi. - ^a^) = (41) 
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The first two terms give — fca(O)^ — — a;^(0) — — which we set to zero 
as an initial condition on ujaix). More precisely if cua is constant it must equal the 
equilibrium value for the given mass. Then we are left with two equations 

fc^(o)(i^^^-r!„.^) = o 

- n^^) + (F^ - Qaf = (42) 

4 Solutions 

In order to look for solutions we set fl^ = 0. The first part of ( ^2| ) gives 

kf^mX = (43) 

If we assume spatial homogeneity then this equation reduces to 

dtf:, = 0. (44) 

Now we want to ask what this result implies about the distribution functions of the 
original degrees of freedom. Using (^2|) we obtain, 

Coj-I- cos^ 6 uj^ sin^ ^ \ / ^ \ 
I 1=0 (45) 

UJ- COS 9 ui- sm ^ j \ Ot 12 / 

In general, this implies that the distribution functions in terms of the original de- 
grees of freedom must also be constant. The only way to obtain non-trivial solutions 
to this equation is to require that the determinant of the matrix of coefficients must 
vanish. This requires 

cos^9sm^eLjlu'i(^-^]=0. (46) 

Since cos 9 and sin 9 can only vanish when the coupling, nip, vanishes (cf (po|)) the 
only possible solution to this equation is for uji to equal uj2. 



5 A Generalizable Analysis 

An explicit rotation of variables is not the approach normally used in solving field 
theories, since it does not readily generalize to arbitrary higher order interactions. 
An alternative method, which leads to an effective theory for one of the systems 
can be found by making a change of variable: 

$2 = (j)2- j {dx')G2{x,x')J2{x') 

= cj)2-mlJ{dx')G2{x,x')<fi{x'). (47) 

The shift involves a non-local Green function, which is assumed to satisfy 

(-□ + ml)G2{x, x') = S{x, x'). (48) 

The boundary conditions on 6*2(0;, x') are not determined by this equation or by 
the Schwinger action principle itself. To determine these boundary conditions we 
use the fact that the interacting fields <i>i and $2 are infiuenced by one another only 
after they are coupled at time t = 0. This physical consideration leads us to choose 
the retarded Green function. 
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Making the shift of variables (^7|) in (|ll|) gives, 

(da;){$t(a;)(-n + m?)^! (x) - j {dx')^l{x)G2{x,x')^i{x') 

+ c^*2{x){-U+ml)c^2{x)] (49) 

We are left with a non-local interaction in terms of the original variables. The 
resulting equations of motion have the matrix form, 

a(*;)=o (50) 

with 

/ 4 \ 

-□ + ml - 



0=\ ^^'"1 ^+^, " (51) 



-□ 



2 / 

To find the particular part of the corresponding Green function we need to invert 
the matrix (D. We obtain, 

detO I, -□ ) 

Using 

detO = (-□ + rn\){-U + mf} - rrip = (-□ + m\){-n + mi) (53) 

it is clear that the Green function has poles at the excitation masses m^ and r7i_. 
This result agrees with the result found in the previous section: the true mass shell 
of the excitations is given by the eigenmasses of the system. 

In order to find the complementary part of the Green function we follow the 
method of the previous section. We define the matrix 

A = ( -ll ^ ) (54) 

The equation that is analogous to ([l^) is then 

$ = A(/). (55) 

We can obtain the diagonalized matrix O given in ( pT] ) by performing the transfor- 
mation on O analogous to (p^: 

d = {A^)-'OA-K (56) 

Using O as defined in (|lj) we reproduce (^). The complementary piece of the 
Green function is obtained from the equation analogous to (p^): 

T={A^)-^TA-K (57) 

To proceed we rewrite (^) in the form 



* " (-□ +ml} ' 
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and impose OT = by writing 

\ f2lS{X2) /22(5(X2) J ^ ' 

where S{x±) indicates any linear combination of i5(x+) and d{x-)- Using ( ]3l| ) and 
(54) in ( |57| ) and comparing with ( |59| ) we obtain, 

/22(5(X2) = f2S{x2) 

4 

/ii<5(x±) = /i-Jlxi) + 7^^-^^/2<5(x2) (60) 

In order to study the time dependence of the distribution functions !F we proceed 
as in the previous section. The (^2(2;) Wightman function obeys the equation of 
motion, 

(-□ +ml)Gi'^\x,x') = 0, (61) 
and its dispersion constraint is the free mass-sheU 

X2 = + ml=0. (62) 
We postulate an ansatz of the form, 

Gi+\x,x') = 2m [ (rfk)e^'=*+^r^^-''^''ii±M^H|g)lf)). (63) 
J 2^22(2:) 

Substituting into (^ij) gives, 

at/22 = (64) 



Using (|62) we get 81/2 — or /2 — constant. 

The $1(2:) Wightman function is more complex; it satisfies the equation of 
motion. 



(-□ +ml)G\~'\x,x')~m* j {dx")G2{x,x")GY\x" ,x') = Q . (65) 

To simplify this expression we operate on the left with (— □ +1112) which gives, 

[(- □ +m^)(- □ +ml) - m2]Gj+^ = (- □ +m\){- □ +m^)Gl+^ = (66) 
We postulate an expression for G'J^'' (x, x') of the form, 

M+)t t\ n ■ [ fj, ^ iki+i dJ^" (1 + /ll(^ll(S),x)) 

G\ '{x,x) = 2m / dk e Jo — . 67 

J 2luii{x) 

Substituting ( |67|) into ( |66| ) we could solve for the dispersion relation, as in the 
previous section. The formalism discussed in this section could be used to study a 
system with a more complicated non-linear interaction, which would lead to non- 
trivial solutions. 
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6 Conclusion 



We have studied a system which consists of two initially independent scalar systems 
that are minimally coupled at some initial time. We have used the machinery 
of Green functions and effective field theory to study the time evolution of the 
statistical distribution functions. This problem is an inherently non-equilibrium 
one. We have discussed how to parameterize such an equilibrating system, and 
how to perform calculations when equations of motion and thus dispersion relations 
change with time. Building on the earlier work in ||^, ^ we have adopted a covariant 
approach to handle time dependent dispersion relations and used the underlying 
conformal structure to show the consistency of working with perturbation theory in 
which dispersion relations are time- independent. We have shown that by formally 
introducing covariant connections, one can rescale a slowly varying non-equilibrium 
theory so that it appears to be an equilibrium one, for the purposes of calculation. 
This approach turns out to be of only formal utility in the example provided, but can 
be generalized to more complex and realistic models. The simple minimally coupled 
system studied here does not reach a conventional kind of macroscopic equilibrium. 
In order to study more conventional equilibration, in which the statistical properties 
are related to the mechanical interactions, we would need to include higher order 
interactions. 
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